Exploiting the co-crystal ligands shape, features and structure-based approaches for identification of SARS-CoV-2 Mpro inhibitors.

SARS-CoV-2 enters the host cell through the ACE2 receptor and replicates its genome using an RNA-Dependent RNA Polymerase (RDRP). The functional RDRP is released from pro-protein pp1ab by the proteolytic activity of Main protease (Mpro) which is encoded within the viral genome. Due to its vital role in proteolysis of viral polyprotein chains, it has become an attractive potential drug target. We employed a hierarchical virtual screening approach to identify small synthetic protease inhibitors. Statistically optimized molecular shape and color-based features (various functional groups) from co-crystal ligands were used to screen different databases through various scoring schemes. Then, the electrostatic complementarity of screened compounds was matched with the most active molecule to further reduce the hit molecules' size. Finally, five hundred eighty-seven molecules were docked in Mpro catalytic binding site, out of which 29 common best hits were selected based on Glide and FRED scores. Five best-fitting compounds in complex with Mpro were subjected to MD simulations to analyze their structural stability and binding affinities with Mpro using MM/GB(PB)SA models. Modeling results suggest that identified hits can act as the lead compounds for designing better active Mpro inhibitors to enhance the chemical space to combat COVID-19.Communicated by Ramaswamy H. Sarma.

Computational investigation of potent inhibitors against SARS-CoV-2 2'-O-methyltransferase (nsp16): Structure-based pharmacophore modeling, molecular docking, molecular dynamics simulations and binding free energy calculations.

The Coronavirus Disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has created unprecedented public health and economic crises around the world. SARS-CoV-2 2'-O-methyltransferase (nsp16) adds a "cap" to viral RNA to maintain the stability of viral RNA, and inhibition of nsp16 activity may reduce viral proliferation, making this protein an attractive drug target. Here, we report the identification of several small molecule inhibitors of nsp16 by virtual screening. First, the nsp16-sinefungin complex (PDB ID: 6WKQ) was selected from the protein data bank. Asp6912, Cys6913, Asp6897 and Asp6928 were determined to be the key amino acids for sinefungin binding in the crystal structure of nsp16-sinefungin complex by molecular dynamics simulation. The complex structures in the stable binding trajectory of nsp16-sinefungin were than clustered through molecular dynamics RMSD analysis. Six clusters were generated, and six representative structures were selected to construct the pharmacophore based on the structure. These six pharmacophores were superimposed on the binding pocket to simplify and pick the common characteristics. The compounds obtained by the pharmacophore screening from Bionet and Chembiv databases were docked into the nsp16 active pocket. The candidate compounds were selected according to the molecular docking score and then screened by MM/GBSA. Finally, four candidate compounds were obtained. Four sets of 150ns molecular dynamics simulations were performed to determine whether candidate compounds could maintain stable interactions with key amino acids. The results of MD and MM/PBSA energy decomposition indicated that C1 and C2 could form a stable complex system with nsp16, and could form strong hydrogen bonds and salt bridges with the key amino acid Asp6897 and Asp6928. This study thus identifies and attempts to validate for the first time the potential inhibitory activities of C1 and C2 against nsp16, allowing the development of potent anti-COVID-19 drugs and unique treatment strategies.